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1 Introduction 

We present a new biomagnetic instrument, designed 
specifically for detection of electrophysiological ac¬ 
tivity within the pregnant maternal abdomen, in¬ 
cluding signals from fetal brain (fMEG), fetal heart 
(fMCG), and uterine smooth muscle. This instru¬ 
ment has been given the acronym “SARA” (SQUID 
Array for Reproductive Assessment). Its principal 
function is to measure the fetal MEG. The magnetic 
fetal auditory evoked response was first detected by 
Blum et al. [1]. Despite improvements in SQUID 
instrumentation, fetal brain signals have remained 
elusive and difficult to observe. Although the fetal 
MEG has been confirmed by several laboratories [2- 
5], results have been disappointing, due to poor sig- 
nal-to-noise ratio. In general, the fMEG measure¬ 
ments are complicated by interference from the fetal 
and maternal MCG, close proximity of the fetal 
heart and head, dependence of the fMEG signal 
maxima upon head position [3], and the fetus motion 
during data acquisition. These factors favour mini¬ 
mizing measurement time, while maximizing the 
chance that there are at least some sensors placed 
close to the fetal signal maxima. 

2 Methods 

2.1 Instrumental design 

The design of an instrument for fMEG involves con¬ 
siderable compromise. Our principal choice was 
between a gantry-based dewar, where a small array 
of sensors is suspended over the mother’s abdomen, 
or a fixed sensor array with substantial coverage. A 
gantry-based system was deemed to have several 
disadvantages. A small array could be designed to fit 
to most areas of the abdominal surface - except for 
the perineal area. However, multiple probe positions 
are needed to map the fetal response. Furthermore, 
in a small signal space, it becomes more difficult to 
resolve the fetal MEG from the fetal MCG. The pos¬ 
sibility of fetal movement and a prolonged meas¬ 
urement session are obvious disadvantages of a 
small array. A larger area gantry-based system pre¬ 
cludes mapping below the pubic symphysis, where 



Figure 1: Cutaway view, showing the front of the 
SARA dewar, reveals its array of 151 primary 
sensors, comprised of 1 st -order axial gradiometers. 

fetal signals may be largest (in late gestation after 
the head is engaged in the pelvic girdle). Lastly, 
with a suspended sensor array, the mother must be 
supine, which is uncomfortable to maintain for ex¬ 
tended periods. 

The most logical alternative to the gantry-based 
system is a stationary, floor-mounted system where 
the mother sits and leans her abdomen against a 
suitably shaped sensing surface. This design is in¬ 
herently safe, the mother is comfortable, and can 
easily mount and dismount the system. Having the 
mother seated on the dewar has the potential disad¬ 
vantage of adding vibrational noise from the mother 
to the sensors. However, we have already demon¬ 
strated a solution to vibrational noise using high-or¬ 
der synthetic gradiometers, with little or no attenua¬ 
tion of the desired signals [6]. 

In the completed design, an array of 151 SQUID 
sensors is organized for coverage of the mother’s 
anterior abdominal surface, from the perineum to the 
top of the uterus (in late gestation). Primary sensor 
flux transformers are axial l st -order gradiometers, 
with 8 cm baselines. The nominal SQUID noise 






Figure 2: A patient is shown seated on SARA. The 
patient support chair slides back to allow the 
patient to easily mount and dismount the instrument. 
The auditory stimulus tubing and light pipe for 
photic stimulus can be seen fastened to the patient’s 
right side by an elastic belt. 

density is 4 fT/VHz. The primary sensor array 
(cutaway in Fig. 1) is curved to fit the pregnant ab¬ 
domen, covering a region of -45 cm high and -33 
cm wide, with an area of -1300 cm 2 and inclined at 
-45 deg. The dewar operates in horizontal orienta¬ 
tion, and has sufficient capacity for one week op¬ 
eration between liquid helium refills. A set of 29 
reference SQUID sensors is incorporated for at¬ 
tenuation of environmental and vibrational noise [7]. 
Flux transformers and SQUIDs are built into the 
liquid helium reservoir, using a proprietary helium- 
vacuum seal. An adjustable padded seat matches the 
height of the perineal dewar extension, so that the 
mother sits upon and leans against the sensor sur¬ 
face, as shown in Figure 2. 

2.2 Signal analysis and noise reduction 

Environmental and vibrational noise are attenuated 
by synthesis of a 3 rd -order gradiometer response for 
the primary sensors [6]. Multiple adaptive methods 
are available for attenuating the maternal MCG and 
fMCG. Accurate cancellation of the maternal MCG 
can be accomplished using multiple regression in the 
Fourier domain. A suitable reference (either ECG or 
MCG) that is free of fMEG signal is required. 
Regression subtraction cannot be used for 
attenuating the fMCG, because any reference chan¬ 
nel containing fMCG is likely to contain fMEG 
components. In order to subtract the fMCG with the 
least damage to the fMEG, nulls in signal space can 
be computed, using the projection operator P = I - 

X(X T X) _1 X T , where X is a matrix of unit vectors 


pointing in the signal-space directions to be nulled. 
Multiple nulls may be required for adequate cancel¬ 
lation of all fMCG components. Lastly, adaptive 
beamforming, and independent components analysis 
(ICA) [8] were evaluated for separating the signal 
components for fMEG, fMCG, maternal MCG, and 
other coherent noise sources. The preceding noise 
reduction methods are normally followed by signal 
averaging. 

Synthetic aperture magnetometry (SAM) was used 
to image the difference in band-limited source 
power, on a voxel-by-voxel basis, throughout the 
abdomen [9]. We used the pseudo-T statistic, based 
upon: 
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where 0 denotes the target voxel position and cur¬ 
rent vectors, and S 2 and v 2 are estimates of the pro¬ 
jected source power and instrumental noise for each 
target voxel. Source power and noise are derived 
from band limited covariance matrices for the active 
and control MEG data segments. The advantage of 
this technique is that a significant source difference 
can be found, even if the brain is not responding 
synchronously with the stimulus. 

2.3 Preliminary measurements 

To date, we have studied the auditory evoked re¬ 
sponse (AEF) to transient tone bursts, and spontane¬ 
ous measurement of uterine smooth muscle during 
the onset of labor. Two types of auditory studies 
have been performed. One paradigm is designed for 
signal averaging of the AEF response to transient 
tone bursts. The other auditory paradigm was de¬ 
signed exclusively for SAM analysis, having an ac¬ 
tive period with a train of frequent tone bursts fol¬ 
lowed by a silent control period. 

The subjects were pregnant women, from 33 to 37 
weeks gestational age, and all gave informed con¬ 
sent for the procedures. Diagnostic ultrasound was 
used to determine the location and orientation of the 
fetus, and especially the fetal head, prior to fMEG 
measurement. Patients were then seated on the 
SARA instrument, leaning forward against the sen¬ 
sor area of the dewar. A reference frame was pro¬ 
vided by localization coils. Auditory stimuli were 
generated by a loudspeaker mounted external to the 
shielded room, and sound was conducted to the pa¬ 
tient’s abdomen by 1 inch ID flexible plastic tubing 
and a soft-rimmed anesthesia mask (Figure 2). 

The averaged AEF studies used 700 Hz or 1 kHz 
tone bursts, with duration of 50 ms and ISI of 2.0 ± 
0.1 seconds. Output was adjusted for a sound pres¬ 
sure level (SPL) of 120 dB, resulting in an estimated 
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Figure 3: Fetal AEF response to 1 kHz tone bursts 
at 33 weeks gestation (overlay of all channels; av¬ 
erage of 310 trials (filtered to DC-30 Hz). 

SPL within the abdomen of between 70 and 80 dB. 
MEG data were recorded in a continuous mode at a 
sample rate of 312.5 Hz, in a band from DC to 100 
Hz, for a total of 350 tone bursts. Signal averaging 
was applied, after projecting out (nulling) the ma¬ 
ternal and fetal MCG components. 

The SAM studies for AEF consisted of alternating 5 
sec trials of 50 msec duration, 1 kHz tone bursts (6 
pulses per second, 120 dB SPL at mouth of tubing), 
and 5 sec trials without stimulus. In total, 50 trials 
were collected (25 during tone stimulus and 25 
during silence). Data were acquired at a sample rate 
of 625 Hz in a DC to 150 Hz bandpass. SAM analy¬ 
sis proceeded first by computing two separate co- 
variance matrices - one for all active and one for all 
control trials. Two volumetric images were gener¬ 
ated from the covariance matrices - a difference im¬ 
age of active and control trials, and an image of the 
sum of active and control trials. 

Signals were collected from smooth uterine muscle 
during early labor at a sample rate of 250 Hz, in a 
DC - 80 Hz bandpass. We collected up to 28 min¬ 
utes of data, divided into epochs of 4 minutes dura¬ 
tion. Data were lowpass filtered (at 1.0 Hz) for 
viewing the low frequency uterine smooth muscle 
activity. For spectral analysis, an FFT was applied to 
segments of uterine activity. 

3 Results 

An averaged fetal AEF was observed in some, but 
not all patients. An observed fetal auditory response 
is shown in Figure 3. A peak signal of nearly 80 fT 
is present in this example, with a latency of 250 ms 
to the first peak. The field map indicated one field 
maximum in the lower abdomen, and the other in 
the perineal area. Ultrasound revealed the fetus to be 
in a normal head down position, with the head en¬ 
gaged in the pelvic girdle. 

A SAM analysis of the same data set as in Fig.3 


(shown in Figure 4), displays a statistically signifi¬ 
cant region of suppression of 4 to 40 Hz activity 
(pseudo-T=3.0). This was consistent with the loca¬ 
tion of the fetal head via ultrasound. A single-state 
SAM image also revealed the position of the fetal 
heart, which also agreed with diagnostic ultrasound. 
Although the signal averaged AEF was not always 
seen, nearly all subjects showed statistically signifi¬ 
cant regions of suppression (active minus control 
power) for the SAM auditory response paradigm. 
Spontaneous measurements of abdominal signals 
during uterine contraction revealed bursts of activity 
including 0.16 Hz (Fig.5). Each burst appears to be 
associated with a contraction, and lasted 
approximately one minute. Although the signal was 
strongest over the left side of the mother’s abdomen, 
it also appeared at the top of the sensor array. These 
signals were not artifacts of respiration. 

Both fetal and maternal MCG were observed with 
high signal-to-noise ratio in all patients. These sig¬ 
nals could be separated using a beamformer. In one 
patient, we observed an elevation in the P-Q seg¬ 
ment of the fetal MCG (Figure 6). 

4 Discussion 

As hypothesized, both the fetal and maternal heart 
signals obscure the fetal averaged AEF in MEG re¬ 
cordings, prior to artifact reduction. The maternal 
MCG signal was easily removed by decorrelation 
subtraction, as it is well separated in signal space. 
Attenuation of the fetal MCG proved more difficult, 
as more than two signal-space nulls were required to 
reduce the fMCG to acceptable levels. Selecting the 
signal-space vector for cancellation required consid- 



Figure 4. SAMpseudo-T statistic image of fetal AEF 
response to tone bursts, at 33 weeks gestation. Each 
voxel represents the source power difference 
between trials of tone bursts and trials of silence, in 
the 4-40 Hz frequency band. The sensor positions 
are shown as a constellation of dots. The bright area 
near the leg cutouts is a region of suppression 
(active minus control), with a pseudo T-value of 3.0. 
This is consistent with the location of the fetal head 
in ultrasound. A marker is also shown for the origin 
of the fMCG signal. 
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Figure 5: Putative uterine smooth muscle signal. 
Sensor plot showing spatial distribution of a 0.16 
Hz spectral line , during early labor. 


erable care, so as not to attenuate the fMEG signal. 
Perhaps our most significant finding is the disparity 
between the averaged evoked response, SAM differ¬ 
ential images, and ICA results. The SAM images 
indicated statistically significant regions of activa¬ 
tion, near the fetal head (from ultrasound examina¬ 
tion). An averaged AEF was visible in only a few of 
these patients. The largest AEF signal was only 80 
fT peak (of 14 patients). Preliminary ICA results 
reveals multiple components that appear to be re¬ 
lated to the fetal AEF. These inconsistencies have 
several potential explanations: 1) poor synchrony of 
the fetal brain’s response to the auditory stimulus as 
a result of poor myelination, or 2) activation of 
multiple pathways. Either mechanism could result in 
an attenuated averaged response, while at the same 
time show significant differences in cortical power 
during the active (stimulated) and control (resting) 
time windows. We hope to resolve these mecha¬ 
nisms, in our pursuit of quantitative fetal MEG. In 
conclusion, both the instrumentation and analysis 
methods have proved successful. 



Figure 6: Interfering signals may be separated from 
the remaining signals by projection operators. 
Note that, with response to DC, a P-Q segment 
elevation can be seen in the fetal MCG. 
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